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Abstract 

A detailed study of the ground state of the Coulomb system (aaee) which corresponds to the 
He^ molecular ion in a magnetic field B = — 4.414 x 10 13 G in parallel configuration (infinitely 
massive a— particles are situated along a magnetic field line) is presented. The variational method 
is employed using a trial function which includes electronic correlation in the form exp ("fr 12) where 
7 is a variational parameter. It is shown that the quantum numbers of the lowest total energy 
state depend on the magnetic field strength. It evolves from the spin-singlet 1 S 5 metastable state 
at < B < 0.85 a.u. to a repulsive spin-triplet 3 T, U state for 0.85 a.u. < B < 1100 a.u. and, finally, 
to a strongly bound spin-triplet 3 II U state at stronger fields B > 1100 a.u. 
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I. INTRODUCTION 



It is well-known that very strong magnetic fields may appear on the surface of white 
dwarfs, B ss 10 6 — 10 9 G and neutron stars, B 10 11 — 10 13 G (for latter case they can 
be even stronger). It seems natural to expect that in a vicinity of the surface atomic and 
molecular systems can occur. For many years this assumption motivates a development of 
atomic-molecular physics in a strong magnetic field since it might lead to understanding of 
the spectra of these compact stars. 

So far, a major attention was paid to a study of one-electron atomic and molecular systems 
(for s review see e.g. [1]). During the last years it was predicted the existence of many exotic, 
strongly bound one-electron molecular system. In turn, two-electron systems had explored 
very little being mostly focused to a study of atomic-type systems H~,He with the only 
exception of the H 2 molecule (see [8], [9] and references therein). Just recently, a first detailed 
study of the molecular ion H% was carried out [2]. It manifested quite surprising evolution 
of the ground state as a magnetic field increases. A goal of the present paper is to make a 
first study of the Coulomb system (aaee) in a magnetic field and establish the existence of 
the molecular ion He% + and its excited states. It is worth mentioning that the molecular 
ion Hel + exists in metastable state in field- free case (see e.g. [3]) being characterized by a 
well-pronounced potential well at a finite internuclear distance. 

Atomic units are used throughout (h—m e —e—l), although energies are expressed in Ry- 
dbergs (Ry). The magnetic field B is given in a.u. with B = 2.35 x 10 9 G. 

II. GENERALITIES 

Let us consider the Coulomb system {aaee) placed in a uniform constant magnetic field. 
We assume that the a— particles are infinitely massive (Born-Oppenheimer approximation 
of zero order). They are situated along the magnetic field line (parallel configuration). The 
Hamiltonian which describes this system when the magnetic field is oriented along the z 
direction, B = (0, 0, B) is written as follows 
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where f)£ = —iVe is the 3-vector of the momentum of the £th-electron, the index k runs over 
the a— particles A and B, r± 2 is the interelectron distance and S = §1 + s 2 is the operator of 
the total spin. At is a vector potential which corresponds to the constant uniform magnetic 
field B chosen in the symmetric gauge, 

At = ^(Bx T t ) = ^(-y t , x e , 0) . (2) 

After substitution of (2) to (1) we arrive at the Hamiltonian in the form 

" - f H + T p ') + ^ + l + *L. + 2S.) , (3) 

i=i v 7 e,K 

where L z = l Zl + l Z2 and S z = s Zl + s Z2 are the z-components of the total angular momentum 
and total spin, respectively. The variable pe = \/ x\ + y\ is the distance from £th electron 
to z— axis. 

The problem under study is characterized by three integrals of motion: (i) the operator of 
the z-component of the total angular momentum (projection of the angular momentum on 
the magnetic field direction) giving rise to the magnetic quantum number m, (ii) the spatial 
parity operator P(fi — > — r\,ri — > — r<i) which has eigenvalues p = ±l(gerade/ungerade) 
(iii) the operator of the z-component of the total spin (projection of the total spin on 
the magnetic field direction) giving rise to the spin quantum number m s . Hence, to any 
eigenstate three explicit quantum numbers can be assigned: the magnetic quantum number 
m, the parity p and the spin quantum number m s . Therefore, the space of eigenstates is 
split into subspaces (sectors) with each of them being characterized by definite values of m, 
p and m s . It is worth noting the Hamiltonian 7i is also invariant with respect to z\ — > —z\ 
and Z2 — > — z 2 (z-parity operator P z ) with quantum numbers <7n = ±1 for positive/negative 
z-parity (this symmetry accounts for the interchange of the nuclei A and B). 

Thus, to classify eigenstates we follow the convention widely accepted in molecular physics 
using the quantum numbers m, S and p. In particular, the notation is 2S+1 M p , where 
2S + 1 is the spin multiplicity and corresponds to 1 for singlet (S = 0) and 3 for triplet 
(S = 1), for the label M Greek letters S,n, A are used that correspond to the states with 
\m\ = 0, 1,2, respectively, and the subscript p (the spatial parity quantum number) will 
take gerade/ungerade(g/w) labels describing positive p = +1 and negative p — — 1 parity 
states. However, there exists a relation between the quantum numbers corresponding to the 
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z-parity (interchange of nuclei A and B) and the spatial parity: 

p=(-l) m a N . 

We restrict our consideration to the states with m = 0,-1, —2 because the ground state of 
a sector with m > always has larger total energy than those with m < 0. 

As a method to explore the problem the variational procedure is used. The recipe of choice 
of trial functions is based on arguments of physical relevance (see e.g. [5]). Eventually, a 
trial function for the state with magnetic quantum number m is chosen in a form 

^) = (1+ <Tc P 12 )(1 + (t jv P A b) 

where a e = 1,-1 stands for spin singlet (S = 0) and triplet states (S = 1), respectively, 
when o"at = 1,-1 stands for nuclear gerade and ungerade states, respectively. The Pyi is the 
operator which interchanges electrons (1 2) and Pab is the operator which interchanges 
the two nuclei A and B. The parameters ai_4, /3i-2 and 7 are variational parameters. If the 
internuclear distance R is taken into account as a variational parameter the trial function 
(4) depends on eight parameters. 

Calculations were performed using the minimization package MINUIT from CERN-LIB. 
Multidimensional integration was carried out using a dynamical partitioning procedure: the 
domain of integration was divided into subdomains following an integrand profile and then 
each subdomain was integrated separately (for details see e.g. [2]). Numerical integration 
was done with a relative accuracy of ~ 1CT 6 — 1CT 7 by use of the adaptive D01FCF routine 
from NAG-LIB. Computations were performed on a dual DELL PC with two Xeon processors 
of 2.8 GHz each (ICN), 54-node FENOMEC and 32-node TOCHITL clusters (UNAM) and 
DUKE dual DELL PC with two Xeon processors of 3.06 GHz each (CINVESTAV). 
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TABLE I: Total Et, double ionization Ej energies (in Ry) and equilibrium distance (in a.u.) of 
the He^ ion in the state Magnetic fields for which Hef^ does not exist (see text) are marked 
by (*). Total energy Ex(He + (lso) + He + (lso)) (in Ry) for zeroing total electron spin projection 
from [2] is shown for a comparison. 
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Q See Ref.[10] 
III. RESULTS 

It is known that in field-free case the He% + molecular ion exists in the state *E S (S = 0, 
o"jv — 1 and m = 0) as a metastable state - it decays to He + + He + (see e.g. [10]). A natural 
goal is to check the existence of bound state of the system (aaee) in parallel configuration 
in presence of the magnetic field ranging B = — 4.414 x 10 13 G. To carry out such a 
study the trial function (4) with a e — 1, ctjv — 1 and m = is used, which depends on 
8 variational parameters. The calculations indicate clearly the existence of a minimum in 
the total energy curve Et{R) for all magnetic fields B > 0. The concrete results for the 
X S 9 state are presented in Table I. This state is a short-lived state for weak magnetic fields 
where its total energy curve lies above the total energy in the dissociation channel to two 
separate helium ions E T (He + (ls ) +He + (ls )), for which the spins of electrons in both He + 
are antiparallel to the magnetic field direction. In general, as the magnetic field strength 
grows, the total energy increases, the system becomes more bound - the double ionization 
energy Ej = 2B — E T increases and more compact (the internuclear equilibrium distance 
decreases) . 

As a next step a detailed study for the 3 S U (5' = 1, ctjv = — l,m = ) state of the Hel + 
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FIG. 1: B = 0: total energies of the l Y, g (S = 0, ground state) and 3 E U (S = 1) states vs. the 
internuclear distance R. 1 Y tg (S = 0) decays to two He + ions, their total energies are shown by 
dashed line. 

molecular ion in parallel configuration in the domain of magnetic fields < B < 4.414 x 
10 13 G is carried out. For this state the variational trial function (4) is characterized by 
<7 e = —1, <7n = — 1 and m — 0, it depends on 8 variational parameters. 

The performed variational calculations do not indicate to any minimum in the total 
energy curve E T (R) of the [pta.ee) system at finite internuclear distances for magnetic fields 
ranging B = — 4.414 x 10 13 G. Thus, the state 3 S M is a repulsive state for all the range of 
studied magnetic fields, < B < 4.414 x 10 13 G. One can state that this system "exists" 
in a form of two Helium ions He + situated at infinitely large distance from each other. 

However, it can be seen that at B pa 0.85 a.u. the crossing occurs of the total energy 
curve at the equilibrium internuclear distance for the 1 S 9 state and the (repulsive) total 
energy curve of the 3 S U state. It means that the He\ + molecular ion exists in the 1 S g as a 
metastable state for < B < 0.85 a.u. However, for > B < 0.85 a.u. the H e\ + ion 
ceases to exist as a compact system characterized by a finite internuclear distance. But it 
does exist in a form of two separate Helium ions He + with each of them in Is state with 
parallel electron spins. As an illustration we present on Fig. 1-3 the plots of total energies of 
the 1 S 9 and 3 S U states for magnetic fields 0, 0.5 a.u. and 10 a.u. 

As a next step the state 3 n n (S = 1, cr N = l,m = —1) of the (aaee) system in parallel 
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FIG. 2: B = 0.5: total energies of the 1 S 9 (S = 0, ground state) and 3 £ u (S = 1) states vs. the 
internuclear distance R. 1 T lg (S = 0) decays to two He + ions, their total energies for different total 
electronic spins are shown by dashed lines. 

configuration is studied in the domain of magnetic fields 2 x 10 9 G < B < 4.414 x 10 13 G. 
For this state the variational trial function (4) with a e = —1, = 1 and m — — 1 depends 
on 8 variational parameters. Obtained results (see Table II) indicate to the existence of a 
clear minimum in the total energy curve Et(R) for magnetic fields B > 10 a.u. With an 
increase of the magnetic field strength the total energy at equilibrium position decreases, the 
system becomes more bound (in this case double ionization energy is equal to Ej = —E T 
and it grows) and more compact (the internuclear equilibrium distance decreases). Also, 
for magnetic fields B > 10 a.u. there exists at least one longitudinal vibrational state 
(see Table II). Table II also shows the total energies of the dissociation channels (aaee) — > 
He + (ls )+He + (2p_i), (aaee) — > He\ + (la g )+e and (aaee) — > He(^(—l)~ l )+a for different 
magnetic fields. 

In Table II the total energy of two Helium ions E T (2) = 2E T (He + (ls )) is shown in the 
case of parallel spins of electrons. This energy also corresponds to the lowest total energy 
of the 3 S U state (see above). The lowest total energies of the 3 I1 U and 3 S U states cross 
at B rs 1100 a.u. At larger magnetic fields the lowest total energy of the system (aaee) 
in the 3 I1 U state is lower than the lowest total energy of the 3 S U state (which correspond 
to the separate ions He + . Hence, the system (aaee) corresponds to the molecular ion 




7 




FIG. 3: B = 10: total energies of the X S S (S = 0) and 3 S U (5 = 1, "ground state") states vs. the 
internuclear distance R. 1 S 9 (5 = 0) decays to two He + ions, their total energies for different total 
electronic spins are shown by dashed lines. 



He\ + at B > 1100 a.u. One can draw a conclusion that the molecular ion Hel + exists at 
B > 1100 a.u. It is a stable system where its ground state is given by the 3 I1 U state. At 
smaller magnetic fields 0.85 a.u. < B < 1100 a.u. the lowest total energy of the system 
(aaee) corresponds to the state 3 S U at the infinite equilibrium distance. It means that the 
system (aaee) does not form a bound state, it does "exist" in a form of two separate Helium 
ions He + . 

If the system (aaee) can be kept externally in the 3 II U state, one can check a sta- 
bility of the system towards dissociation channels. A comparison of the total energy 
of (aaee) in the 3 U U state and the total energy of the (dominant) channel of decay 
(aaee) — > He + (ls ) + He + (2p_ 1 ) indicates that the (aaee) system in the 3 U U state be- 
comes stable with respect to the possible channels of decay for B > 100 a.u. The dissociation 
energy increases monotonously with a magnetic field growth. 

To complete a study of the existence of the Hel + ion we carried out calculations of 
the energy for the state 3 A 9 (S — 1, gn — 1, m — —2) of the system (aaee) in parallel 
configuration in the domain of magnetic fields 1 a.u. < B < 4.414 x 10 13 G. For this state 
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TABLE II: The state 3 II M : total Et, lowest longitudinal vibrational E v n, energies (in Ry) and 
equilibrium distance R eq (in a.u.) of the He^ ion as a function of the magnetic field. Magnetic 
fields for which Hel + does not exist (see text) are marked by (*). Total energy of two Helium 
ions He + (lso) + He + (2p_i) with electron spins antiparallel to magnetic field (from [2]) is denoted 
as Et(1), also the total energy of two Helium ions both in ls-state, He + (lso) each of them with 
electron spin antiparallel to magnetic field from [2] is denoted as Et(2) for comparison as well 
as total energy of the system He^il&g) + e from [2]. Total energy ET{He(*)) for He atom in 
3(— l) _1 )-state was taken from [7]. 
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TABLE III: The state 3 A g : total energy Et in Ry and equilibrium distance R eq in a.u. of the 
Hel + ion. Magnetic fields for which He^ does not exist (see text) are marked by (*) . 
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the variational trial function (4) is used with a e = —1, <jn — 1 and m = —2, it depends on 
8 variational parameters. 

There exists a clear minimum in the total energy surface Et{R) of (aaee) for magnetic 
fields B > 10 a.u. It is found that with an increase of the magnetic field strength the total 
energies decrease, the system becomes more bound (double ionization energies increase, 
Ei = —Et, for triplet states) and more compact (the internuclear equilibrium distance 
decreases) (see Table III). At B w 7000 a.u. the crossing between the total energy curves 
of the 3 S U and 3 A g state occurs. The 3 A 9 state becomes the lowest-lying excited state for 
B > 7000 a.u. until the Schwinger limit. 
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IV. CONCLUSION 



The existence of the Hel + molecular ion and its low-lying excited states in a magnetic 
field ranging from zero to 4.414 x 10 13 G in parallel configuration are studied in the Born- 
Oppenheimer approximation using the variational method. It is found that the quantum 
numbers of the state of lowest total energy depend on the magnetic field strength. The 
ground state evolves from the spin-singlet l Y, g state which is a metastable state for weak 
magnetic fields B < 0.85 a.u. to the unbound (repulsive) spin-triplet 3 £ u for intermediate 
fields 0.85 a.u. < B < 1100 a.u. and eventually to a strongly bound spin-triplet 3 II M state for 
strong magnetic fields B > 1100 a.u. It manifests the existence of the stable He\ + molecular 
ion at B > 1100 a.u. 

In the domain 1100 a.u. < B < 7000 a.u. where the state 3 n u becomes the ground 
state, the lowest-lying excited state is the repulsive 3 S M state. The next excited state is 
the "bound" 3 A g state (having a clear minimum in the total energy curve) lies above the 
3 S M state. One can state that the He\ + molecular ion does not have excited states in a 
traditional sense in this domain of magnetic fields. At larger magnetic fields B > 7000 a.u. 
up to the Schwinger limit the lowest-lying excited state is 3 A 9 , while the next excited state 
is the repulsive 3 S M state. 
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